Voter model dynamics in complex networks: Role of dimensionality, disorder and 

degree distribution. 
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We analyze the ordering dynamics of the voter model in different classes of complex networks. We 
observe that whether the voter dynamics orders the system depends on the effective dimensionality 
of the interaction networks. We also find that when there is no ordering in the system, the average 
survival time of metastable states in finite networks decreases with network disorder and degree 
heterogeneity. The existence of hubs in the network modifies the linear system size scaling law 
of the survival time. The size of an ordered domain is sensitive to the network disorder and the 
average connectivity, decreasing with both; however it seems not to depend on network size and 
degree heterogeneity. 

PACS numbers: 64.60.Cn,89.75.-k,87.23.Ge 



I. INTRODUCTION 

Equilibrium order-disorder phase transitions, as well 
as nonequilibrium transitions and the kinetics of these 
transitions have been widely study by spin Ising- 
type models in different lattices 0- Given the recent 
widespread interest on complex networks 0. [H 
the effect of the network topology on the ordering pro- 
cesses described by these models has also been considered 
0, H H El EH In particular models of opinion forma- 
tion, or with similar social motivations, have been dis- 
cussed when interactions are defined through a complex 
network [H [H Hi [H [H [H . 

A paradigmatic and simple model where a systematic 
study of network topology effects can be addressed is 
the voter model [l9|. for which analytical and well es- 
tablished results exist in regular lattices [2(J, |21|. The 
dynamics of ordering processes for the voter model in reg- 
ular lattices |22| is known to depend on dimensionality, 
with metastable disordered states prevailing for d > 2. 
In this paper we address the general question of the role 
of network topology in determining if the systems or- 
ders or not, and on the dynamics of the ordering process. 
Specifically, analyzing the voter model in several different 
networks, we consider the role of the effective dimension- 
ality of the network, of the degree distribution and of the 
level of disorder present in the network. 

The paper is organized as follows. In section 2 we 
shortly review the basics as well as recent results on 
the voter model. Section 3 considers the voter model 
scale free (SF) networks 0] of different effective 



(SF 



dimensionality 34], showing that voter dynamics can or- 
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der the system in spite of a SF degree distribution. In 
section 4 we consider the role of network disorder by in- 
troducing a disorder parameter that leads from a struc- 
tured (effectively one-dimensional) SF (SSF) network |23| 
to a random SF (RSF) network through a small world [24j 
SF (SWSF) network. The role of the degree distribution 
is discussed comparing the results on the SSF, RSF and 
SWSF networks with networks with an equivalent disor- 
der but without a power law degree distribution. Some 
general conclusions are given in Section 5. 



II. VOTER MODEL 

The voter model is defined by a set of "voters" 
with two opinions or spins Oi = ±1 located at the nodes 
of a network. The elementary dynamical step consists 
in randomly choosing one node (asynchronous update) 
and assigning to it the opinion, or spin value, of one 
of its nearest neighbors, also chosen at random. In a 
general network two spins are nearest neighbors if they 
are connected by a direct link. Therefore, the probability 
that a spin changes is given by 




(1) 



where fcj is the degree of node i, that is the number of 
its nearest neighbors, and Vi is the neighborhood of node 
i, that is the set of nearest neighboring nodes of node 
i. In the asynchronous update used here, one time step 
corresponds to updating a number of nodes equal to the 
system size, so that each node is, on the average, updated 
once. In our work we choose initial random configura- 
tions with the same proportion of spins +1 and — 1. 

The dynamical rule implemented here corresponds to 
a node-update. An alternative dynamics is given by a 
link-update rule in which the elementary dynamical step 
consists in randomly choosing a pair of nearest neighbor 
spins, i.e. a link, and randomly assigning to both near- 
est neighbor spins the same value if they have different 
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FIG. 1: Interface density evolution for individual realizations, 
for system size N = 10000 and average connectivity (k) = 8. 
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FIG. 2: Plateau height for SF networks with different average 
connectivity. Network size N = 10000. Inset: Plateau height 
dependence on network size. Average connectivity (k) — 6. 
Data is averaged over 1000 realizations. 



values, and leaving them unchanged otherwise. These 
two updating rules are equivalent in a regular lattice, 
but they are different in a complex network in which dif- 
ferent nodes have different number of nearest neighbors 
pfij . In particular, both rules conserve the ensemble av- 
erage magnetization in a regular lattice, while in a com- 
plex network this is only a conserved quantity for link- 
update dynamics. Node- update dynamics conserves an 
average magnetization weighted by the degree of the node 
. We restrict ourselves in this paper to the stan- 
dard node-update for better comparison with the grow- 
ing literature on the voter model in complex networks 

[USE! Ill 

The voter model dynamics has two absorbing states, 
corresponding to situations in which all the spins have 
converged to the <7j = 1 or to the Oi = — 1 states. The 
ordering dynamics towards one of these attractors in a 
one-dimensional lattice is equivalent to the one of the 
zero temperature kinetic Ising model with Glauber dy- 
namics. In more general situations, as in regular lattice 
of higher dimension or in a complex network, the ordering 
dynamics is still a zero-temperature dynamics driven by 
interfacial noise, with no role played by surface tension. 
A comparison of the voter model and the zero tempera- 
ture Ising Glauber dynamics in complex networks has 
been recently reported [30j . A standard order parameter 
to measure the ordering process in the voter model dy- 
namics [22,|22| is the average interface density p, defined 
as the density of links connecting sites with different spin 
value: 

/ n \ N 

0- EE^P /!>• m 

\j=i jeVi ) »=i 

In a disordered configuration with randomly dis- 
tributed spins p ~ 1/2, while when p takes a small value 
it indicates the presence of large spatial domains in which 
each spin is surrounded by nearest neighbor spins with 



the same value. For a completely ordered system, that 
is, for any of the two absorbing states, p — 0. Starting 
from a random initial condition, the time evolution of p 
describes the kinetics of the ordering process. In regular 
lattices of dimensionality d < 2 the system orders. This 
means that, in the limit of large systems, there is a coars- 
ening process with unbounded growth of spatial domains 
of one of the absorbing states. The asymptotic regime of 
approach to the ordered state is characterized in d = 1 by 
a power law (p) ~ while for the critical dimension 
d = 2 a logarithmic decay is found (p) ~ (lni) -1 [2^ . 
Here the average (•) is an ensemble average. 

In regular lattices with d > 2 j2fj, as well as in small 
world networks , it is known that the voter dynamics 
does not order the system in the thermodynamic limit 
of large systems. After an initial transient, the system 
falls in these cases in a metastable partially ordered state 
where coarsening processes have stopped: spatial do- 
mains of a given attractor, on the average, do not grow. 
In the initial transient of a given realization of the pro- 
cess, p initially decreases, indicating a partial ordering of 
the system. After this initial transient p fluctuates ran- 
domly around an average plateau value £. This quantity 
gives a measure of the partial order of the metastable 
state since I = gives an estimate of the average lin- 
ear size of an ordered domain in that state. In a finite 
system the metastable state has a finite lifetime: a fi- 
nite size fluctuation takes the system from the metastable 
state to one of the two ordered absorbing states. In this 
process the fluctuation orders the system and p changes 
from its metastable plateau value to p — 0. Considering 
an ensemble of realizations, the ordering of each of them 
typically happens randomly with a constant rate. This is 
reflected in an exponential decay of the ensemble average 
interface density 

(p) oc e~~ , (3) 
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FIG. 3: Interface density evolution in SF networks of different 
sizes (increasing from left to right: N = 1000, 2000, 5000, 
10000, 20000, and 50000). Data is averaged over at least 
1000 realizations for (k) = 6. 




FIG. 4: Interface density evolution in one dimensional sys- 
tems with and without scale-free topology. For reference, the 
dotted line is a power law with exponent —^.Average over 
1000 realizations, N = 10000 and (k) = 8. 



where r is the survival time of the partially ordered 
metastable state. Note then that the average plateau 
value £ has to be calculated at each time, averaging only 
over the realization of the ensemble that have not yet 
decayed to p = 0. 

The survival time r, for a regular lattice in d = 3 
pcj and also for a small world network (2^ , is known to 
scale linearly with the system size N, r ~ N, so that 
the system does not order in the thermodynamic limit. 
More recentl y t he same scaling has been found for ran- 
dom graphs |2^, l3pj w hile a scaling r ~ jv 88 has been 
numerically found |25ll3til | for the voter model in the scale 
free Barabasi- Albert network . This scaling is com- 
patible with the analytical result r ~ N/ In N reported in 
Ref. |29j. Other analytical results for uncorrelated net- 
works with arbitrary power law degree distribution are 
also reported in Ref. |2||. We note that a conceptually 
different, but related quantity, is the time that a finite 
system takes to reach an absorbing state when coarsen- 
ing processes are at work. This time n is known to scale 
as T\ ~ N 2 for a regular d = 1 lattice and t\ ~ iVlniV 
for a regular d = 2 lattice. 

In the next sections we discuss the time evolution of p 
and the characteristic properties of the plateau value £ 
and survival time r for asynchronous node-update voter 
dynamics in a variety of different complex networks. 

III. DIMENSIONALITY AND ORDERING: 
VOTER MODEL IN SCALE-FREE NETWORKS 

One of simplest models that displays a scale-free de- 
gree distributions is the well known Barabasi- Albert net- 
work [32]. In this model, the degree distribution follows 
a power law with an exponent 7 = 3, the path length 
grows logarithmically with the system size [3j while the 
clustering coefficient decreases with system size j^- It 



has been shown that critical phenomena on this class of 
networks are well reproduced by mean field calculations 
valid for random networks |33| . Thus we will consider in 
the remainder the Barabasi- Albert networks as a repre- 
sentative example of a random scale-free (RSF) network. 
Results for the voter model in the BA network are shown 
in Figs. 11131 The qualitative behavior that we observe 
is the same than the one described above for regular lat- 
tices of d > 2 or also observed in a small world network 
[27j : The system does not order but reaches a metastable 
partially ordered state. The interface density p for dif- 
ferent individual realizations of the dynamics is shown in 
Fig. ^ In this figure we see examples of how finite size 
fluctuations take the system from the metastable state 
with a finite plateau value of p to the absorbing state 
with p = 0. The level of ordering in this finite lifetime 
metastable state can be quantified by the plateau level 
£ shown in Fig. [21 We find that the level of ordering 
decreases significantly with the average connectivity of 
the network, a result consistent with the idea that total 
ordering is more easily achieved for effective lower dimen- 
sionality. On the other hand the level of ordering is not 
seen to be sensitive to the system size, for large enough 
sizes. 

The survival time r can be calculated from the ensem- 
ble average interface density (p) as indicated in Eq. 
The time dependence of (p) for systems of different size 
(Fig. |3) shows an exponential decrease for which the re- 
sult mentioned above r ~ ]\[0.88 can ^ e obtained [2o| . 
We note that the value r is found to be independent of 
the mean connectivity of the network and that a linear 
scaling r ~ N is obtained if a link update dynamics is 
used j25j. 

The fact that the presence of hubs in the BA network 
is not an efficient mechanism to order the system might 
be counterintuitive, in the same way that the presence 
of long range links in a small world network is also not 



efficient to iead to an ordered state. However, in both 
cases the effective dimensionality of the network is infin- 
ity and the result is in agreement with what is known for 
regular lattices with d > 2. A natural question is then 
the relevance of the degree distribution versus the effec- 
tive dimensionality in the ordering dynamics. To address 
this question we have chosen to study the voter model 
dynamics in the Structured Scale Free (SSF) network in- 
troduced in Ref. |2^|. The SSF networks are a nonran- 
dom network with a power law degree distribution with 
exponent 7 = 3 but with an effective dimension d = 1 

0. 

Our results for the time dependence of the average in- 
terface density in the SSF network are shown in Fig. 01 
For comparison the results for a regular d = 1 network 
are also included. For both networks we observe that 
the system orders with the average interface density de- 
creasing with a power law with characteristic exponent 
1/2 



(p) (t) ~ t 



(4) 



The only noticeable difference is that the SSF network 
has a larger number of interfaces at any moment, but the 
ordering process follows the same power law. Addition- 
ally we find that for a finite systems the time t\ to reach 
the absorbing state scales as t\ ~ ^V 2 , as it also happens 
for the regular d = 1 network: 

The network is completely ordered when the last in- 
terface disappears. At this point, the density is simply 
(N (fc)) _1 , where ./V (k) is the total number of links in the 
network. Since the interface density decreases (p) ~ t~?, 
then the time to order t± is given by 



-1/2 



(5) 



leading to n ~ N 2 . 

Therefore we conclude that the effective dimensionality 
of the network is the important ingredient in determin- 
ing the ordering process that results from a voter model 
dynamics, while the fact that the system orders or falls 
in a metastable state is not sensitive to the degree distri- 
bution. 




FIG. 5: Mean interface density evolution for SWSF networks 
with different disorder parameter p (increasing from circles 
to stars: p = 0.0004, 0.001, 0.002, 0.004, 0.01, 0.02, 0.04, 
and 0.1), and comparison with the BA network (+). For 
reference, the dotted line is a power law with exponent — i. 
Data averaged over 100 realization for p < 0.01 and over 1000 
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FIG. 6: Survival times for a SWSF network with different dis- 
order parameter p. For comparison results for SW networks 
are also included. Average over 1000 realizations, TV = 10000 
and (k) = 8. 



IV. ROLE OF NETWORK DISORDER AND 
DEGREE HETEROGENEITY 

Once we have identified in the previous section the 
crucial role of dimensionality we now address the role of 
network disorder and degree heterogeneity in quantita- 
tive aspects of the voter model dynamics. We do that by 
considering a collection of complex networks in which the 
system falls into partially disordered metastable states, 
except for the the regular one-dimensional lattice and 
SSF networks in which the system shows genuine order- 
ing dynamics: 

1. Structured Scale Free (SSF) network as defined in 
the previous section. 



2. Small-World Scale-Free (SWSF) network. This is 
defined by rewiring with probability p the links of a 
SSF network. In order to conserve the degree dis- 
tribution of the unperturbed (p = 0) networks, a 
randomly chosen link connecting nodes i, j is per- 
muted with that connecting nodes k, I) |3q . 

3. Random Scale Free (RSF) network: Defined as the 
limit p = 1 of the SWSF network. The RSF net- 
work shares most important characteristics with 
the BA network. 

By changing the parameter p from p = (SSF) 
to p = 1 (RSF) we can analyze how increasing 
levels of disorder affect the voter model dynam- 
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FIG. 7: Plateau heights for SWSF networks with different 
disorder parameter p. For reference, dotted line is a power 
law £ ~ p while the solid line £ ~ p 76 . For comparison 
results for SW networks are also included. Average over 100 
realization for p < 0.01 and over 1000 realizations for p > 
0.01, N = 10000 and (k) = 8. Inset: Plateau heights for 
SWSF networks of different system size N. Average over 1000 
realizations, with p — 0.01 and (k) = 8. 



ics while keeping a scale free degree distribution. 
On the other hand, the consequences of the de- 
gree heterogeneity characteristic of SF networks 
can be analyzed comparing the voter model dynam- 
ics on these networks with networks with the same 
level of disorder and a non-SF degree distribution. 
These other networks are constructed introducing 
the same disorder parameter p, but starting from a 
regular d = 1 network. Namely we consider: 

4. Regular d = 1 network that can be compared with 
a SSF network. 

5. Small World (SW) network defined introducing the 
rewiring parameter p in the regular network as in 
the prescription by Watts and Strogatz |24|. The 
SW network can be compared with the SWSF net- 
work. 

6. Random network (RN) corresponding to the limit 
p = 1 of the SW network. 

Likewise, one can consider a random network with 
an exponential (EN) degree distribution. The EN 
network is constructed as in the BA prescription 
but with random instead of preferential attachment 
of the new nodes. These two random networks, RN 
and EN, can be compared with the RSF network. 



A. Role of disorder 

Figure[S]shows the evolution of the mean interface den- 
sity for SWSF networks with different values of the dis- 
order parameter p. It shows how varying p one smoothly 




FIG. 8: Interface density evolution for SFSW and SW net- 
works with the same level of disorder p = 0.01. Average over 
1000 realizations, N = 10000, and (k) = 8. 



interpolates between the results for the SSF network and 
those for a RSF network. In general, increasing net- 
work randomness by increasing p the system approaches 
the behavior in a BA network, making it to to fall in a 
metastable state of higher disorder, but with finite size 
fluctuations causing faster ordering. This trend is quanti- 
tatively shown in Fig.[S]and Fig.[7|where the the survival 
time t and plateau level £ for SWSF networks are plot- 
ted as a function of the disorder parameter p. We observe 
that r and the size of the ordered domains I — de- 
crease with p but without following any clear power law. 
As a general conclusion, when extrapolating to p = 1, we 
find that t S wsf > trsf and 1 S wsf > Irsf- 

The role of increasing disorder in the network can also 
be analyzed in networks without a scale free degree distri- 
bution by considering SW networks with different values 
of the rewiring parameter p. The survival time r and 
plateau level £ for SW networks are also plotted in Fig. El 
and Fig. [7\ We observe that the effect of disorder is 
qualitatively the same for SW than for SWSF networks 
[35| . Extrapolating the results in Fig. Eland Fig. to 
p = 1 where the SW network becomes a RN we find that 
tsw > trn and Isw > Irn- 



B. Role of degree distribution 

To address the question of the role of the degree dis- 
tribution of the network in the voter model dynamics 
we compare the evolution in networks with a scale free 
degree distribution with the evolution in equivalent net- 
works but with a degree distribution involving a single 
scale. A first comparison was already made between the 
dynamics in a regular d — 1 network and the SSF net- 
work (Fig. 2J). This is included for reference in Fig. |S| 
where we compare the evolution of the mean interface 
density in a SWSF network with the evolution in a SW 
network with the same level of disorder. We observe for 
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FIG. 9: Mean interface evolution for SWSF networks of dif- 
ferent system size N (increasing from left to right: N = 1000, 
2000, 5000, 10000, 20000, 50000). Average over 1000 realiza- 
tions, with p — 0.01 and (k) — 8. 



SWSF a similar plateau value (similar but slightly more 
disordered state) at any time before the exponential de- 
cay of (p) which is faster for the SWSF than for the 
SW networks. Finite size fluctuations that order the sys- 
tem seem to be more efficient when hubs are present, 
causing complete ordering more often, and therefore a 
faster exponential decay of (p). These claims are made 
quantitative in Fig. [S] and Fig. [7| where it is shown that 
t~sw > tswsf and Isw — Iswsf- In addition, extrapo- 
lating to the limit p = 1 we have that trn > trsf and 
Irn — Irsf 

It is also interesting to compare the dependence with 
system size of the voter model dynamics in SW |27| and 
SFSW networks: The time dependence of the mean in- 
terface density for a SWSF network with an intermediate 
fixed value of p is shown in (Fig. EJ). The qualitative be- 
havior is the same than the one found for SW networks. 
However, the survival times shown in Fig. 1101 deviates 
consistently from the linear power law r ~ N found for 
SW networks jS^. This deviation might possibly have 
the same origin than the deviation from the linear power 
law observed for BA networks, that is the lack of con- 
servation of magnetization in the node-update dynamics 
of the voter model in a complex network |25| . This non- 
conservation becomes much more important in a SWSF 
network than in a SW network because of the high degree 
heterogeneity. On the other hand we note that the an- 
alytical results for survival times in Ref. 29| apply only 
to uncorrelated networks and therefore do not help us in 
understanding our numerical result for SWSF networks. 
We also mention that the plateau level £ for SWSF net- 
works does not show important dependence with system 
size (see inset of Fig. 0) 

The role of degree heterogeneity can be further clari- 
fied considering the limit of random networks p = 1 where 
the SW network becomes a RN and the SWSF network 
becomes a RSF network essentially equivalent to the BA 




FIG. 10: Survival times for SWSF networks of different sys- 
tem size N. Average over 1000 realizations, with p — 0.01 
and (k) = 8. 



network. The evolution for the mean interface density 
for different random networks is shown in Fig. 1111 We 
find again that when there are hubs (large degree het- 
erogeneity) there is a faster exponential decay of (p), so 
that ordering is faster in BA networks than in RN or EN, 
while the plateau level or level of order in that state does 
not seem to be sensitive to the degree distribution . This 
coincides with the extrapolation to p — 1 of the data in 
Fig. Eland Fig. [7| which indicates that trn > trsf and 
Irn — Irsf- Our results for the system size dependence 
of the survival times and plateau levels for RN and EN 
networks are shown in Fig. The size of the ordered 
domains I = is again found not to be sensitive to 
system size. The survival times for RN and EN networks 
follow a linear scaling r ~ N in agreement with the pre- 
diction in Ref. 29]. We recall that, as discussed earlier, 
in random networks with scale free distribution such as 
the BA network a different scaling is found (r ~ JV ' 88 ) 
[HHii compatible with the prediction r ~ N/\nN p9j. 



V. CONCLUSIONS 

We have analyzed how the ordering dynamics of the 
voter model is affected by the topology of the network 
that defines the interaction among the nodes. First we 
have shown that the voter model dynamics orders the sys- 
tem in a SSF network [2^, which is a scale- free network 
with an effective dimension d = 1. This result, together 
with the known result that in regular lattices the voter 
model orders in d < 2, suggests that the effective dimen- 
sion of the underlying network is a relevant parameter 
to determine whether the voter model orders, but not 
its degree distribution. The relevance of the effective di- 
mensionality of different scale-free networkshas also been 
observed in other dynamical processes 0, HE • 
In the SSF network the density of interfaces in the voter 
model decreases as (p) ~ t~ x /' 2 in the same as in the 
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FIG. 11: Mean interface density evolution for RN and EN of 
different sizes N. BA network of size 10000 is also shown for 
comparison. Average over 1000 realizations, and (k) — 8. 
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FIG. 12: Survival times for RN and EN networks of different 
sizes N. Average over 1000 realizations, with (k) = 8. Inset: 
Plateau heights for the same networks. 



metastable disordered states so that the survival time to 
reach an ordered state in finite networks is smaller 

TSWSF > trsf, t S w > Trs ■ 

Likewise, the average size of ordered domains in these 
metastable states decreases with increasing disorder 



SWSF 



> I 



RSF, 



hw > I 



RS 



Third, the degree heterogeneity also facilitates reach- 
ing an absorbing ordered configuration in finite networks 
by decreasing the survival time: finite size fluctuations 
ordering the system are more efficient when there are 
hubsin the network, so that 

T SW > T SWSF, Tun > TjiSF ■ 

The presence of hubs also invalidates the scaling law for 
the survival time r ~ N fouund in SW and RN. However 
we didn't appreciate differences in the average size of 
ordered domains depending on degree heterogeneity 



sw 



Iswsf, Irn — Irsf 



In summary, we find for the different classes of net- 
works considered in this work that 

TSW > T RN > T R SF 

Isw — Iswsf > Irn — Irsf ■ 

In general our results illustrate how different features 
(dimensionality, order, degree heterogeneity) of complex 
networks modify key aspects of a simple stochastic dy- 
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one-dimensional regular lattices. 
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